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I. INTRODUCTION 

The considerable interest in the properties of heavy- 
fermion antiferromagnets is due to their unconven¬ 
tional superconductivity, quantum phase transitions, 
pronounced competition between a tendency to magnetic 
ordering and Kondo fluctuations, and coexistence of su¬ 
perconductivity and antiferromagnetism. The quantum 
phase transitions are initiated by an external or chem¬ 
ical pressure and are accompanied by the variations in 
the ground state structure, which leads to the change in 
the characteristics of materials. In particular, in heavy- 
fermion CeCue-ajAua; and YbRh 2 Si 2 metals, a passage 
through the quantum critical point is accompanied by 
destruction of the long-range antiferromagnetic (AFM) 
order varying the control parameters, specifically, dopant 
concentration x and magnetic held [1]. Under pres¬ 
sure, the phase diagrams of the Ce-based compounds, 
including CePd 2 Si 2 , Celna [2], CeRhIns [3, 4], and 
CePt 2 ln 7 [5], contain a superconductivity dome in the 
vicinity of the expected quantum critical point. 

The nature of magnetic ordering is one of the most im¬ 
portant problems of physics of heavy-fermion systems. If 
we assume that the long-range AFM order is initiated by 
the RKKI indirect exchange interaction and Kondo huc- 
tuations tend to destruct the magnetic ordering, then at 
such a scenario of the competition [6] at the quantum 
critical point, simultaneously with the destruction (oc¬ 
currence) of antiferromagnetism, the Kondo regime can 
be established (suppressed) [7, 8] and the transition from 
localized to delocalized electrons can occur [9]. 

According to the modern concepts, 4f electrons in the 
AFM phase of cerium compounds are quasi-localized and 
form the coherent heavy-fermion state. This was con¬ 
firmed by the experimental data, which show that even in 
the CeRhIns AFM phase the effective and cyclotron elec¬ 
tron masses are larger than the free electron mass [10, 11]. 
According to study [12], the mixed-valence regime can be 
implemented in this compound. Therefore, it is reason¬ 
able to investigate the formation of magnetic ordering 
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using the periodic Anderson model (PAM) in the regime 
when the localized bare level is close to the Fermi level. 
Such an approach allows describing the strong renormal¬ 
ization of electron mass and the Fermi surface topology 
variation at a quantum critical point without using the 
Kondo breakdown scenario [13-15]. 

The occurrence of the AFM phase in the PAM 
was demonstrated, e.g., in the Hatree-Fock approxima¬ 
tion [16] and using the slave-boson technique [17, 18]. It 
should be noted that in the considered approaches Neel 
temperature Tn takes small values only around the quan¬ 
tum phase transitions (see also [19]). It means that even 
minor variations in the external factors, e.g., pressure, 
can induce the quantum phase transition from the AFM 
to paramagnetic (PM) phase. However, in many heavy- 
fermion compounds with Tn of no more than few Kelvin 
degrees, antiferromagnetism is sufficiently stable against 
pressure variations. 

The alternative approach to finding a magnetic insta¬ 
bility point is the calculation of dynamic magnetic sus¬ 
ceptibility. In the PAM, the dynamic magnetic suscepti¬ 
bility was calculated using equations of motion for irre¬ 
ducible Green’s functions [20] and in the random phase 
approximation [21]. In study [22], it was proposed to 
calculate the dynamic magnetic susceptibility using the 
perturbation theory for the hybridization interaction. It 
was demonstrated that in the mixed-valence regime the 
effective interaction caused by hybridization between lo¬ 
calized and itinerant electrons suppresses any magnetic 
fluctuations. In the limit [/ —>■ oo, where U is the pa¬ 
rameter of on-site Coulomb interaction, the method for 
calculating the dynamic magnetic susceptibility was de¬ 
veloped on the basis of a diagram technique for Hubbard 
operators within the Hubbard model and t—J model [23]. 
In study [24], this method was applied to determine the 
dynamic magnetic susceptibility in the PAM paramag¬ 
netic phase. 

Heavy-fermion compounds, e.g., Ce„TmIn 3 „+ 2 m [25], 
have a quasi-two-dimensional (quasi-2D) structure 
schematically illustrated in Fig. 1. As is known, the Neel 
temperature of a quasi-2D Heisenberg antiferromagnet 
is approximately determined as Tn = ! {}-n{J/K) + c), 

where J is the parameter of exchange between the nearest 
ions in the (xy) plane, K is the value of exchange interac- 
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tion between the nearest neighbors along the z axis, and 
c is the constant depending on the lattice type [26] . This 
formula indicates a decrease in the transition temper¬ 
ature relative to the isotropic case. In study [27], using 
the neutron spectroscopy data and Heisenberg model, the 
parameters of exchange between Ce ions in the CeRhIns 
compound were found to be J = 0.74 meV and K = Q.\ 
meV. 



FIG. 1. Schematic of the quasi-two-dimensional structure of 
Cells compounds. J and K are the constants of magnetic 
exchange between Ce ions. 


It should be noted that in the most interesting case of 
strong electron correlations for the electron states of rare- 
earth ions and the one-electron excitation energy close 
to the Fermi level (the mixed-valence case), it is conve¬ 
nient to divide the hybridization processes into high- and 
low-energy ones [28]. The high-energy processes are the 
transitions at which, due to the strong correlations, the 
energy of a system changes by a value much larger than 
the hybridization interaction parameter. This great en¬ 
ergy difference makes it possible to take into account the 
above-mentioned hybridization mixing using the opera¬ 
tor perturbation theory in the atomic representation and 
to obtain the effective Hamiltonian describing, in particu¬ 
lar, the exchange coupling between quasi-localized states 
of rare-earth ions. The parameter of this interaction is 
determined as J ^ where V is the intensity of 

hybridization between localized and itinerant electrons. 
The rest low-energy contributions determine the proper¬ 
ties of the mixed-valence regime. Note that in the Cell5 
compounds (e.g., CeRhIns), the role of localized electrons 
is played by Ce 4f electrons and the collective states are 
formed mainly by In p electrons [29]. 

In this work, using the effective PAM, which explic¬ 
itly takes into account the exchange interaction between 
4f electrons, we obtain a pressure dependence of the 
Neel temperature for quasi-2D cerium intermetallic com¬ 
pounds. This dependence is not only qualitatively consis¬ 
tent with the data reported in [30] , but also describes well 
the experimental results. The pressure dependence of 
the Neel temperature consists of two portions. The first 
portion is characterized by a linear decrease in the Neel 


temperature with increasing pressure. Such a behavior is 
qualitatively reproduced with disregard of hybridization. 
The second portion shows a sharp break of antiferromag¬ 
netism. We show that this dependence is formed only 
with regard to the low-energy hybridization processes. 
In this case, two channels affecting the magnetic order¬ 
ing occur: the exchange interaction of 4f electrons tends 
to establish the AFM ordering, whereas the low-energy 
hybridization of f and p electrons can both promote and 
suppress antiferromagnetism. In this study, we estimate 
partial contributions of these microscopic mechanisms to 
the experimental effective parameters of exchange inter¬ 
actions in the CeRhIns compound [27]. 


II. MODEL AND METHOD 


The Hamiltonian of the effective PAM, which takes 
into account the exchange interaction between 4f elec¬ 
trons for a simple cubic lattice, is 


HeG — ^kc\fyCk(T + + 

k(7 raa 

+ J E + 

m^l ^ ' 


( 1 ) 


The first term in the Hamiltonian describes a subsys¬ 
tem of itinerant electrons (In p electrons) in the quasi¬ 
momentum space with bare energy counted from 
chemical potential /i. Localized 4f electrons correspond¬ 
ing to Wannier cell m are described in the atomic 
representation using the Hubbard operators = 

jn; TO)(m; sj, where jn;m) is the atomic states without f 
electron (jO; m)) and with one f electron (jcr; m)) and with 
different spin moment projections a. From the doublet 
in the crystal field, the f level with j = hl2 is taken into 
account. For the bare energy of 4f electron ^lcj, the self- 
consistent held is introduced. The exchange interaction 
is specihed by parameter Jmi, Sm is the quasi-spin vec¬ 
tor operator of f electron, and Nm is the operator of the 
number of localized electrons on site m. The fourth term 
in the Hamiltonian determines the low-energy hybridiza¬ 
tion processes between localized and itinerant electrons 
with intensity 14; A is the number of sites in the lattice. 

To describe the magnetic properties, we apply a dia¬ 
gram technique for the Hubbard operators [31, 32], which 
will be used for calculating the Matsubara Green’s func¬ 
tion in the atomic representation [33] 


D± (mr; ui't') 


(T,Xt4(r)xi1;(r')^(/3)) . (2) 

\ / 0,c 


The time-dependent Hubbard operators are taken in the 
interaction representation; is the time-ordering opera¬ 
tor. Index 0 indicates that averaging is made with regard 
to the Hamiltonian that describes noninteracting itiner¬ 
ant and localized electrons. Only connected diagrams are 
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FIG. 2. General view of quasiparticle spectrum in the an¬ 
tiferromagnetic phase of the periodic Anderson model with 
regard to hybridization and hoppings only in the [xy) plane. 

taken into account. The scattering matrix has the form 
S{(3) = Trexp^—f^ Hint{T)dT^, where interaction op¬ 
erator Hint involves the operators of hybridization and 
exchange interactions and P = 1/T is the inverse tem¬ 
perature. 

As is known [34], the Fourier image of Green’s func¬ 
tion can be presented in the form [q) = Gj_ (g) P (g), 
where g = are the even Matsubara fre¬ 

quencies and P (g) is the force operator. Then, the Dyson 
equation for the function (g) is 

Gi (g) = Gf (g) + cf (g) E (g) (g), (3) 

Gf (g) = G(°) (zw^) + G(°) {lujm) P (g) J,-Gf (g) ,(4) 


exchange due to the residual part of the hybridization 
interaction and corrections 5P (g) renormalize the bare 
vertex. 

The quasi-2D structure of the Cell5 cerium com¬ 
pounds is schematically shown in Fig. 1. Parameter J 
of the exchange between the nearest Ce ions in the {xy) 
plane is significantly larger than the analogous parame¬ 
ter K along the 2 axis. In study [27], for CeRhIns the 
weak exchange between next-to-nearest neighbors along 
the 2 axis was also taken into account, which allowed de¬ 
scribing the incommensurate magnetic structure at am¬ 
bient pressure. In this study, we limit the considera¬ 
tion to the account for the exchange interaction only be¬ 
tween the nearest neighbors, since under pressure the 
CeRhIns AFM structure becomes commensurate [35]. 
Then, the long-range exchange parameters can be ig¬ 
nored. It should be noted that the low-energy mixing 
of f and p electrons, as well as and hoppings of itinerant 
In electrons, are limited only by the {xy) plane. 

To describe antiferromagnetism, it is convenient to 
pass to the two-sublattice representation. Then, the bare 
energy of localized electrons in the F sublattice, the mag¬ 
netization of which is codirectional to the 2 axis, is de¬ 
termined as pFcr = Eq — y — (2J -|- K){n\,l2 -|- r\aR)- 
Here, nn is the f-electron concentration and R = (-S’!), 
ya = 1(—1) at a =t (J.). The bare energy of f electrons 
in the G sublattice is pccj = pFs where d denotes the 
opposite direction of a. In the description of a quasi-2D 
antiferromagnet, all the quantities in Dyson equations 
(3, 4) should be replaced by matrices. Then, the matrix 
Green’s function is determined as 

qAiBi qAiB2\ 
qA^B^ QA2B2 y ^ ^ ^ 


Gx = 


/qFF QFG\ 


ygGF qGGj 


Gi^ = 


where G(°) {iuJm) = {itOm — 2 / 1 ) ^ is the bare Green’s 
function with the self-consistent exchange field h, 
is the Fourier image of the exchange integral, and all 
the corrections related to hybridization are contained in 
mass operator E (g) and force operator P (g). For con¬ 
venience, below we explicitly distinguish the bare vertex 
and corrections related to the hybridization interaction 
P (g) = (5'^) + SP (g) in the force operator. 

The solution of Eq. (3) is given by 

( \ _ _ G^°^ {iUJm) _ , . 

l-[E(g) + J,^P(g)]G(0)(za;^)- 

In the limit R —>■ 0 the solution acquires a simple form 
corresponding to the Tyablikov approximation for the 
Heisenberg model. In the steady-state case, this func¬ 
tion can have the features related to the formation of 
ferromagnetic or AFM ordering at the instability point. 
The exchange interaction between 4f electrons, which 
is induced by high-energy hybridization processes, leads 
to antiferromagnetism and the low-energy hybridization 
processes affect the initial AFM state. In this case, mass 
operator E (g) indicates the occurrence of a new effective 


where A, B = F, G. Notations Fn and Gn {n = 1,2) in¬ 
dicate that this Green’s function is built on the operators 
belonging to the F or G sublattice and corresponding to 
the n plane of the unit cell presented in Fig. 1. The 
matrices for the mass and force operator components are 
written in a similar way. Taking into account the ex¬ 
change parameters denoted on Fig. 1 an interaction ma¬ 
trix is given by: 



/ = 


(J<r kA 
\K<r J<r) ’ 


( 7 ) 


where O is a null matrix. 

Figures 3 and 4 show the general view of diagrams 
for arbitrary components AnBn of the matrix mass and 
force operators. Note that the expressions are indepen¬ 
dent of number n of the plane in the quasi-2D unit cell, 
so below this index will be omitted. In the figures, solid 
lines with two arrows [> and ► indicate propagators in 
the Hubbard-I approximation for localized electrons with 
spin moment projections f £md respectively, with re¬ 
gard to hybridization. The bare Green’s functions for f 
electrons are indicated by the solid line with one arrow 
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t> (►). The solid line with two arrows indicates any 
of the four bold propagators for itinerant electrons in the 
two-sublattice representation. Symbols o and • indicate 
the Hubbard vertex factors Foa = + X^')^ for the 

corresponding electron spin directions. The summation 
is made over the internal momenta p = (p, iwn), where 
are the odd frequencies. The total number of diagrams 
for the matrix mass and force operators is 64. 



FIG. 3. Diagrams for component AB of the mass operator. 


q 




FIG. 4. Diagrams for component AB of the force operator. 


The analytical expression for mass operator component 
FF obtained from the diagrams in Fig. 3 is 




p(T 


\Vp + Wpf GTip) + [Vp - Wpf Gf (p)+ 

+ (W^^-VS) (G“'5(p) + Gf“(p))]. 


( 8 ) 


The functions G^^ (p) and G^ (p) {A, B = F, G\ v, p = 
a, /3) are the propagators for localized and itinerant 
electrons in the Hubbard-I approximation in the two- 
sublattice description, respectively, and Vp and Wp are 
the Fourier images of hybridization integrals inside a sub¬ 
lattice and between sublattices, respectively. 

Using the matrix Dyson equations, it is easy to obtaine 
the Green’s function the poles of which determine 

the spin-wave excitation spectrum for a quasi-2D anti- 
ferromagnet with regard to hybridization processes be¬ 
tween localized and itinerant electrons. These poles are 
described well by approximate analytical expressions: 


^1.2 q — Woi,2 {Y) + q, (9) 

where a;oi ,2 {q) are the bare branches of the spectrum 
of a quasi-2D Heisenberg antiferromagnet (without hy¬ 


bridization) and the corrections 

SiUi ^ = diPF + diCG ^{KiFG + KiOf) - — X 

Wo* yq) 

X {(4J -I- 2K) R [diCG — diFF i (KiPG + Kicp)] + 

+ {J^+ K^)R [diFG — d*GF ± {KiFF + KiGG)]} (10) 

describe the effect of hybridization interaction on the 
magnon spectrum. We introduced the designations 

d,AB = m /2 + (a;o* {q)) /2, (11) 

K,ab = (wo* m /2, (12) 

where the mass and force operator components are cal¬ 
culated after the analytical continuation and F = G, 
G = F. 

The antiferromagnetic order parameter is defined as 
R = ni,l2 — (X|l), where 

Using the obtained self-consistent equation, the Neel 
temperature in the limit i? —> 0 is determined. 


III. RESULTS AND DISCUSSIONS 

A. Spin-wave spectrum of quasi-two-dimensional 
heavy-fermion antiferromagnets 

It follows from the obtained expressions for the spin- 
wave spectrum that the Goldstone theorem about the ex¬ 
istence of gapless magnon excitations in the AFM phase 
at t/ —> oo with regard to hybridization of f and p elec¬ 
trons is valid. The occurrence of two branches in the 
magnon spectrum is caused by different characters of ro¬ 
tation of spin moments on the nearest sites in the {xy) 
plane and along the z axis with regard to the quasi-two- 
dimensionality in the antiferromagnet (inphase and an¬ 
tiphase, respectively). The branch indicated by the 
superscript in formula (10) is Goldstone. For the branch 
ui 2 q, the excitations are separated by an energy gap. 

We consider the case when the Fermi level lies close 
to Eq and crosses the weak-dispersion region in Fig. 2, 
where the general structure of the quasi-particle spec¬ 
trum in the PAM is shown. Then, heavy fermions can 
be formed in the magnetically ordered phase the mass of 
which exceeds the mass of free electrons. We suggest that 
this state is implemented in CeRhIns near atmospheric 
pressure. 

The spin-wave spectrum for the principle direction of 
the AFM Brillouin zone, which corresponds to the f- 
electron concentration ul ~ 0.7, is shown in Fig. 5. 
In the figure, qa is the wave vector component and 
a = X, y, z, one of the unit cell parameters. Dashed 
lines indicate magnon energies Wpi {q) for the quasi-2D 
structure with disregard of the hybridization between p 
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FIG. 5. Spin-wave spectrum of the quasi-two-dimensional 
structure with regard to hybridization between localized and 
itinerant electrons (solid lines) and bare spectrum of local¬ 
ized electrons (dashed lines) for the (111) principal direction 
of the antiferromagnetic Brillouin zone. Concentration of 4f 
electrons is ul ~ 0.7 (see text for details). 

and f electrons. Solid lines are the branches of the 
spin-wave spectrum with regard to the hybridization in¬ 
teraction. It can be seen that the low-energy hybridiza¬ 
tion processes lead to a significant increase in the spin- 
wave stiffness n for the Goldstone mode tuiq = Kq (at 
small q values) and magnon energy. 

The interaction parameters for ul ~ 0.7 were chosen 
in the form V = 0.3|ti|, J = 0.004|ti|, K = J/10, where 
ti is the matrix element of hoppings of itinerant electrons 
between the nearest sites {ti < 0), Eq = 1.5ti is the en¬ 
ergy of the f-level, and Ug = 1.2 is the total electron con¬ 
centration. Comparison of the model law of dispersion 
of itinerant electrons and the dispersion dependencies 
for In p electrons in CeRhIns obtained using ab initio 
calculations yields |ti| ~ 0.1 — 0.3 eV [29]. The chosen 
parameters correspond to the Fermi excitation spectrum 
presented in Fig. 6. 

It should be noted that in the regime when the lo¬ 
calized electron subsystem is almost completely filled 
(riL ~ 1) the magnon spectrum almost does not change 
with regard to the hybridization interaction. 

Modification of the spin-wave spectrum is explained by 
the occurrence of the additional effective exchange inter¬ 
action due to low-energy hybridization processes. Indeed, 
the comparison of Dyson equation (3), which contains the 
mass operator, and equation (4), which takes into ac¬ 
count the exchange interaction between 4f electrons and 
vertex factors, shows that the mass operator components 
work as the effective exchange interaction. The exchange 
parameters of the effective interaction between different 
lattice sites can be estimated as 

( 14 ) 

Q 

E (15) 


FIG. 6. Fermi spectrum for the f-electron concentration 
ul « 0.7. Filling level p is indicated by the dashed line. 

where sites / and /' belong to the F sublattice and sites g, 
to the G sublattice. It should be noted that the effective 
interaction occurs only between ions lying in the (xy) 
plane, since the low-energy hybridization processes are 
limited by this plane. Thus, the exchange interaction 
along the z axis with parameter K determines the quasi- 
2D character of the systems under study and exchange 
parameter J in the [xy) plane is renormalized due to 
the hybridization interaction. In this case, the exchange 
between the next-to-nearest neighbors is also formed by 
the expense of the effective interaction in the [xy) plane. 

The effect of hybridization between p and f electrons 
on the characteristics of the AFM phase, including the 
spin-wave stiffness, AFM order parameter, and Neel tem¬ 
perature, is determined by the signs of parameters Af^g 
and Ifj'- At Af^g > 0 and Ifj' < 0, the hybridization 
processes promote the AFM ordering and the spin-wave 
spectrum acquires the form presented in Fig. 5. If the 
magnon energy decreases with regard to hybridization, 
as it happens near the quantum critical point (see the 
next paragraph), then the AFM exchange between sub¬ 
lattices in the [xy) plane weakens, since Af^g < 0, and 
the exchange inside the sublattice satisfies the inequal¬ 
ity Ifj' > 0. Thus, it is important that the introduced 
effective interaction parameters depend on the localized 
level position, density of states, and temperature. 

B. Pressure dependence of the Neel temperature 

The low-energy hybridization processes not only lead 
to the quantitative variation in the parameters of the 
AFM phase, but also qualitatively change the behavior of 
these parameters. Let us consider the dependence of Neel 
temperature Tn on pressure P in heavy-fermion quasi-2D 
Ge-based antiferromagnets, such as GeRhIns. According 
to the experimental data, the Neel temperature linearly 
decreases with increasing pressure in a fairly wide pres¬ 
sure range [30]. At the critical pressure, the Neel tem¬ 
perature turns to zero and the long-range AFM order 
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is destroyed. Note that the Neel temperature in these 
materials is no higher than few Kelvin degrees. 

It is assumed that the pressure growth leads to an in¬ 
crease in energy Eq of a 4f electron on the positively 
charged Ce ion due to enhancing Coulomb interaction 
with the negatively charged environment. Since this in¬ 
teraction, including that between sites, is the strongest 
in these systems, the effect of the growth of Eg prevails 
over the growth of hybridization intensity and hoppings 
with increasing pressure. 

Figure 7 shows the dependence of the Neel temperature 
on the bare energy of 4f electron (pressure). Dots indicate 
the dependence with regard to the low-energy hybridiza¬ 
tion processes between f and p electrons. This curve sep¬ 
arates the regions of implementation of the AFM phase 
(dashed area) and the PM phase. The solid line shows 
the dependence 7 n(P) with disregard of the hybridiza¬ 
tion interaction. It can be seen that in this case the 
Neel temperature linearly decreases with increasing pres¬ 
sure. Such a behavior is related to a linear decrease in 
4f-electron concentration. The renormalized curve also 
contains a linear portion; however, the Neel tempera¬ 
ture grows due to hybridization. It is more important 
that the account for the hybridization interaction leads to 
destruction of antiferromagnetism with increasing pres¬ 
sure and the dependence T'n(T’) becomes consistent with 
the experiment. The quantitative consistency of the re¬ 
sults with the data for CeRhIns is reached if we take 
|ti| ~ 0.14 eV. This estimation is adequate for heavy- 
fermion systems. It should be noted that at the quan¬ 
tum phase transition point in pressure from the AFM to 
PM phase the Fermi surface broadens and the effective 
electron mass strongly grows [15], as was experimentally 
observed in [10]. 

At Eq = —l.bjti], the spin-wave and Fermi spectra 
are presented in Figs. 5 and 6, respectively. The local¬ 
ized electron concentration is nn ~ 0.7. It can be seen 
in Fig. 7 that this point corresponds to the pressure 
similar to atmospheric. Using formula (14), we estimate 
the effective parameter of exchange between the nearest 
ions caused only by the low-energy hybridization pro¬ 
cesses: Ai = 0.0036jtij, which corresponds to the AFM 
exchange. When estimating this parameter, we excluded 
the coefficient that takes into account the vertex factor. 
Then, the total value of exchange interaction between 
the nearest ions in the {xy) plane is written in the form 
Jeff = J -b Ai, where J = 0.004jtij is the parameter of 
the initial exchange interaction in the PAM, which is in¬ 
duced by the high-energy hybridization processes in the 
limit U —^ oo. Therefore, the Neel temperature increases. 
Note that the effective exchange parameters JeS agrees 
well with the parameter estimated in [27]. 

The Neel temperature can be estimated using a simple 
formula similar to the formula for a quasi-2D Heisenberg 
antiferromagnet [26]: 


10 '' 
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FIG. 7. Dependence of the Neel temperature on the energy 
of 4f electrons (pressure) with regard to hybridization (dots) 
and with disregard of it (solid line). The shaded area indicates 
the antiferromagnetic (AFM) phase and the white area, the 
paramagnetic (PM) phase. 
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where c = 3.15 and K is the parameter of exchange be¬ 
tween the nearest ions along the 2 axis. Thus, the main 
effects related to the low-energy hybridization interac¬ 
tion in the linear portion of the dependence TnIP) can 
be approximately reduced to renormalization of the pa¬ 
rameter of exchange between the nearest Ce ions in the 
(xy) plane. As the pressure is increased, np decreases and 
parameter Ai decreases as well. However, at approaching 
the critical pressure, the effective exchange between more 
distant ions, which was not taken into account in (16), 
becomes important. Near the critical pressure, the frus¬ 
trations arise, which are caused by the competition be¬ 
tween the AFM and ferromagnetic exchange. As a result, 
the long-range order disappears. 

The described effects induced by hybridization of itin¬ 
erant and localized electrons in heavy-fermion antiferro- 
magnets will take place also in 3D compounds, e.g., in 
Celna. However, in the description of these compounds, 
it is impossible to limit the consideration to the account 
for hybridization and hoppings only in the (xy) plane. 
Another interesting feature of cerium intermetallic com¬ 
pounds is the coexistence of superconductivity and an¬ 
tiferromagnetism near the quantum critical point. In 
the proposed model, the exchange interaction between 4f 
electrons can induce the Cooper instability [36]. Then, 
the formation of superconductivity near the quantum 
critical point can be unrelated to quantum fluctuations 
and be explained by the fact that antiferromagnetism 
suppresses Cooper pairing. However, these problems lie 
beyond this study. 


IV. CONCLUSIONS 


In 


tiijTtJ 0ff 

ln(Jeff/Ar) -b c’ 


(16) 


Using the periodic Anderson model, we investigated 
the interference of two microscopic mechanisms of the for- 
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mation of exchange interaction between f electrons in the 
quasi-2D heavy-fermion cerium intermetallic compounds. 

The first mechanism is implemented at a large value 
of the intraatomic Coulomb repulsion and is caused by 
the high-energy hybridization between itinerant p elec¬ 
trons and f electrons of rare-earth ions. The intensity 
of exchange coupling determined by this mechanism is 
independent of temperature, concentration of itinerant 
carriers, and position of the chemical potential level. 

The situation is qualitatively different for the second 
mechanism initiated by the low-energy hybridization be¬ 
tween the above-mentioned electron groups. The contri¬ 
bution of these processes in the resulting exchange cou¬ 
pling between Ce ions significantly depends on the Fermi 
level position and density of states of itinerant electrons. 
This conclusion follows from the analysis of the behavior 


of magnetization of the antiferromagnetic sublattice ob¬ 
tained with regard to the contributions of these processes 
and calculated using a diagram technique for Hubbard 
operators. It was demonstrated that the account for the 
second mechanism plays a decisive role in the satisfac¬ 
tory description of the experimental data obtained for 
the quasi-2D cerium systems, e.g. CeRhIns. 
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